Enhancement of Noise-induced Escape through the Existence of a Chaotic Saddle 
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We study the noise-induced escape process in a prototype dissipative nonequilibrium system, the 
Ikeda map. In the presence of a chaotic saddle embedded in the basin of attraction of the metastable 
state, we find the novel phenomenon of a strong enhancement of noise-induced escape. This result 
is established by employing the theory of quasipotentials. Our finding is of general validity and 
should be experimentally observable. 
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Since the seminal treatment of the noise-induced es- 
cape problem by Kramers 1], major progress has been 
made by Onsager and Machlup. They realized that the 
escape process consists of large fluctuations, which are 
very rare, and that the trajectory peaks sharply around 
some optimal (most probable) escape path 2]- Thus, 
despite the stochastic nature of the escape process, the 
escape path is of almost deterministic nature, as other 
paths than the most probable one have an exponentially 
smaller probability. That theory was derived for a small 
noise level 5 — > 0. A review on noise-induced escape in 
equilibrium systems and the most important recent ad- 
vancements is given in j^. 

In the last decade it has been realized, that systems that 
are not in thermal equilibrium or are lacking the prop- 
erty of detailed balance, can give rise to a large variety of 
interesting phenomena in the noise-induced escape prob- 
lem. Only recently experiments on this problem have 
been conducted, using Josephson junctions, electronic 
circuits, lasers, and an electron in a Penning trap 
Some of the most interesting novel theoretical findings 
include a pre-exponential factor of the Kramers rate [5j, 
a symmetry breaking bifurcation of the optimal escape 
path J3| and a distribution of the escape paths origi- 
nating from a cusp point singularity Furthermore, 
the very intriguing phenomenon of saddle-point avoid- 
ance has been discovered 0- For a fluctuating barrier 
the effect of resonant activation has been theoretically 
predicted l9| and experimentally confirmed Also a 
stepwise growth of the escape rate for short time scales 
has been found Recently, an oscillation of the es- 

cape rate in dependence on the friction for a multiwell 
potential was demonstrated 12] . For periodically driven 
systems a number of interesting results has been obtained 
as welL like a resonantly decrease in the activation en- 
ergy [13, a log arithmic susceptibility of the fluctuation 
probability [lj| , time oscillations of escape rate s 1 151 and 
enhancement of escape due to transient chaos |l6||. 
Here we report on a new mechanism of lowering the re- 
quired energy for noise-induced escape (enhancement of 
escape), thus a reduction of the mean first passage time. 
This happens, if a chaotic saddle is embedded in the open 



neighborhood of the basin of a metastable state. Then 
the escape trajectory does not only pass through a sin- 
gle unstable periodic orbit By contrast, it can pass 
through the chaotic saddle, i. e. a geometrically strange, 
invariant, non-attracting set (which is made up of an infi- 
nite number of unstable periodic orbits). The trajectory 
can jump between points of the chaotic saddle with no 
additional activation energy required. The overall low- 
ering of the activation threshold is due to the fact, that 
the escape process consists now of three subsequent steps: 
Firstly, the trajectory jumps on one orbit on the chaotic 
saddle. Secondly, it switches on the chaotic saddle, with- 
out need of input energy, to select the orbit which allows 
the easiest escape, and thirdly it fluctuates from that or- 
bit to the saddle point on the basin boundary. By this 
mechanism, the chaotic saddle is transformed into a dy- 
namically relevant quantity, whereas in noisefree systems 
it is only important for transient behavior. In this way 
the chaotic saddle acts as a 'shortcut'. 
Since noise-induced escape has previously been studied 
using dissipative maps 18], which allow analysis in a 
straitforward way, we demonstrate our findings for the 
Ikeda map 0]. This is an idealized model of a laser 
pulse i n an optical cavity. With complex variables it has 
the form 



Zn+l 



bzn exp 



irj 



1 + \z„ 



(1) 



where Zn — Xn + iyn is related to the amplitude and 
phase of the nth laser pulse exiting the cavity. The pa- 
rameter a is the laser input amplitude and corresponds to 
the forcing of the system. The damping (1 — h) accounts 
for the reflection properties of mirrors in the cavity and 
measures the dissipation. The empty cavity detuning is 
given by k and the detuning due to a nonlinear dielectric 
medium by 77. The Ikeda map gives rise to rich dynami- 
cal behavior, exhibiting for some parameters even highly 
multistable behavior (23 . 

We fix the parameters at a = 0.85, b — 0.9 and k — 0.4 
and vary only 77 in the range 2.6 < rj < 12. For the noise- 
less system two stable states are present. One fixed point 
(state A, in Fig. ^ undergoes a period doubling see- 
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nario and becomes a chaotic attractor at w 5.5 Another 
fixed point (state B, □ in Fig. ^ remains a fixed point 
over the whole parameter range considered. The noise- 
induced escape from state A is investigated. The basin 
boundary separating these two stable states is a smooth 
curve which is build by the stable manifold of the saddle 
point C (* in Fig. separating the two stable states. 
With the above mentioned features of the system, there 
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FIG. 1: Grey dots represent the basin of attraction for the 
fixed point marked with a for rj = 4.1. The other fixed 
point is also depicted (□). The chaotic saddle is shown with 
black dots and the saddle point on the basin boundary is 
marked by ^. 



tials can be derived through a minimization procedure of 
the action of escape trajectories from a Hamilton- Jacobi 
equation [2^ . The action to be minimized has the form 



N-l 



SN[{zi)o 



<i<N\ 



(2) 



1=0 



for the map Zn+i = f{zn) + where a is the standard 
deviation of the additive, Gaussian, white noise term 
With appropriate boundary conditions the infimum of 
this action with respect to N and i along a path is the 
quasipotential <i>. The mean first exit time is then given 
in analogy to Kramer's law: 



(t) ^ exp 



(3) 



whith A<I> defined as the minimal quasipotential differ- 
ence 

A$ := inf {$(y) - $(a) : a G A, y e dG}, (4) 

where A is the attractor and dG is the basin bound- 
ary. In Fig. 121 the quasipotential is shown for 77 = 4.1. 
There is a single peak corresponding to the fixed point 
A and a plateau region of a practically constant quasipo- 
tential, which refiects the chaotic saddle. Employing 



are no unusual effects expected in the noise-induced es- 
cape problem. However, for a critical value of r^c = 3.5686 
an additional period 3 solution close to the fixed point 
(state A) emerges, with a fractal basin boundary between 
these two solutions and a chaotic saddle embedded in 
this fractal basin boundary. It is important to note, that 
the basin boundary between the two fixed points remains 
smooth over the whole parameter range considered here. 
Increasing 77 further, the stable period 3 solution dissa- 
pears in a boundary crisis, yet a chaotic saddle is still 
present beyond the boundary crisis, completely embed- 
ded in the open neighborhood of the basin of the sta- 
ble solution, as can be seen in Fig. ^ for 77 = 4.1. A 
chaotic saddle is a geometrically strange, invariant, non- 
attracting set. It is computed using the PIM-triple algo- 
rithm |2l|. We stress that it is the chaotic saddle, that 
has a remarkable effect on the average escape time from 
the stable fixed point. 

To treat the problem of noise-induced escape we now em- 
ploy the theory of quasipotentials, which gives rigorous 
results on the infiuence of noise on the invariant density 
and the mean first passage time. Quasipotentials have 
been introduced in the mathematical literature for time- 
continuous systems in |22l| and for discrete time ones in 
[23j |. For systems of physical interest, they were first pro- 
posed in 12 51 an d extended to systems with coexisting at- 
tractors in [2g. Discrete systems with strange invariant 
sets were for the first time treated in ^27]. Quasipoten- 




FIG. 2: Quasipotential 4>(a;, y) for the Ikeda map with rj — A.l 
on a 300 x 300 grid. The single peak corresponds to the fixed 
point. Also, an extended plateau at — log <l>(a;, y) ~ 5.0 is 
visible, caused by the chaotic saddle. 

the quasipotential for the noise-induced escape problem, 
the minimum value of y) on the basin boundary has 
to be determined. This is exactly the minimum escape 
energy A^{x,y), since the quasipotential at the stable 
solution (fixed point, periodic orbit or chaotic attractor) 
is zero. The point on the basin boundary, where this 
happens, is generally a saddle point of the system. 



3 



To quantify the escape process with the quasipotential, 
we plot for various values of r] the corresponding mini- 
mal escape energy A^{x,y) in Fig. |31 To elucidate the 
role of the chaotic saddle as the origin of an enhancement 
of noise-induced escape, we also include in the plot the 
value of the height of the plateau in the quasipotential. 
In the framework of quasipotentials, the difference in 
height of the escape energy and the saddle plateau corre- 
sponds to the distance between the basin boundary and 
the chaotic saddle, whereas the height of the saddle is 
related to the distance between the attractor and the 
saddle. 

The mechanism of the escape process is closely connected 
to the existence of an embedded chaotic saddle. It con- 
sists of two steps, namely, a noise-induced fluctuation 
from the attractor (state A) to the chaotic saddle, and 
then from the chaotic saddle to the fixed point (state 
B) via the saddle point on the boundary. The escape 
can also be incomplete, as the trajectory may fall back 
from the chaotic saddle to the attractor. In a success- 
ful escape, the chaotic saddle acts as a 'shortcut', as its 
presence lowers the overall escape energy. This behavior 
seems to be especially pronounced if the chaotic saddle is 
closer to the basin boundary than to the attractor (com- 
pare the region 3.6 < rj < 4.5 of Fig. Let us note 
that for rj = 5.5 it is not clear, if there exists a chaotic 
saddle, which is the case for all other values of 77 > 3.5686 
we have tested. The PIM-triple method, as well as the 
quasipotential, yield for 77 = 5.5 no conclusive result, as 
a chaotic saddle may exist very close to the chaotic at- 
tractor and numerically it is very difficult to distinguish 
between the two. 
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FIG. 3: Minimal escape energy (0) and height of the saddle 
plateau (A) in logarithmic scale versus rj. 



It is important to quantify the influence of the relative 
sizes of the basins of attraction, since it is increasing with 
increasing rj and we are here only interested in the change 
of activation energy caused by the chaotic saddle. The 
distance between the attractor and the saddle point on 



the boundary is usually proportional to the relative size 
of the basin. Both quantities are expected to play a role 
in the stability of the metastable state located in the 
basin, although we are not aware of any theoretical work 
dealing with this relation directly. To compensate for the 
change of escape energy caused by the increase in size of 
the basin of attraction, in Fig. ^ the escape energy is 
divided by the ratio of the size of the basin of attraction 
of state A to the overall area (A -I- B) for 3 different sec- 
tions of the phase space. The sections of the phase space 
decrease from top to bottom. The combination of the 
two quantities, potential height and basin size, yields a 
pronounced minimum at 77 « 5.0 for all 3 curves, thus 
confirming the essential role in lowering the escape en- 
ergy played by the chaotic saddle. The most probable 
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FIG. 4: Average quasipotential height (= escape energy) 
divided by the ratios of the sizes of the basins of attrac- 
tion (# initial conditions A / (A -I- B)). The curves corre- 
spond to a smaller frame of reference from top to bottom, 
with the values x G [-10.0, 10.0], y G [-10.0,10.0] (marked 
with □), X e [-5.0,5.0],?/ e [-7.0,7.0] (marked with A), 
X G [—3.0, 2.0], y G [—3.5,4.0] (marked with ^), respectively. 

escape path 0, H^j for rj — 5.0 is shown, together with 
the chaotic saddle, in Fig. \S\ For this parameter value, 
there is a stable period 4 solution. As can be seen, the 
trajectory jumps at first directly on points of the chaotic 
saddle, moves secondly along points of the chaotic saddle 
for some iterations, until it is thirdly transported close 
to the basin boundary to the saddle point. Since the 
first step (from fixed point to saddle) and the last step 
(from saddle to the basin boundary) are minimal in this 
case, the enhancement is maximal. Other values of the 
Ikeda map, where no chaotic saddle is present, have been 
investigated as well. For these parameter values the ef- 
fect could not be found, and the graphs corresponding 
to Fig. 2]have a strictly monotonic shape. This demon- 
strates that the existence of the chaotic saddle is of cru- 
cial importance for the occurrence of the enhancement of 
noise-induced escape. 

Moreover, the stability of a fixed point is determined by 
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FIG. 5: Escape path and chaotic saddle for r\ — 5.0. The 
saddle point on the basin boundary is shown as H<. 



its eigenvalues. The eigenvalues are found to be A = 0.9 
for the whole range, where it exists. Consequently, the 
linear approximation is of no relevance to the noise- 
induced escape problem, as its range of validity is much 
smaller than the region for the escape, which is the whole 
open set of the basin of attraction shown in Fig. 
To conclude, we have demonstrated the effect of en- 
hancement of noise-induced escape through the existence 
of a chaotic saddle in the open neighborhood of the 
metastable state for the Ikeda map as a parameter is var- 
ied. Employing the theory of quasipotentials, it was pos- 
sible to understand this lowering of the escape threshold. 
We stress that the reported mechanism of the lowering 
of the escape energy is of qualitatively different nature 
from a recently found effect, where also an enhancement 
of noise- induced escape through transient motion (typical 
for chaotic saddles) has been found In this scenario, 
a nonadiabatically, periodically driven system exhibits a 
facilitation of noise-induced interwell transitions. This 
occurs, because the basin boundary becomes fractal and 
the distance between the two states is effectively reduced. 
In the mechanism reported here we always have a smooth 
basin boundary between the two states and the chaotic 
saddle is embedded in the basin of one state, not in the 
basin boundary between the states. The analysis of the 
exact escape path on the chaotic saddle, in contrast to 
the case where the trajectory leaves via a single periodic 
orbit Ji] will be the presented in a much broader hash- 
ion in a future publication l^^. The reported new phe- 
nomenon is of general relevance for many physical and 
chemical problems. It is predicted to occur in a variety 
of systems, and should experimentally be observable. 
We acknowledge A. Hamm and D. Luchinsky for valuable 
discussions and A. Hamm also for the help in program- 
ming. This work was supported by the DFG and INTAS. 
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